Generation of 'Super-Ponderomotive' Electrons due to a non- Wakefield interaction 
between a Laser Pulse and a Longitudinal Electric Field 
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It is shown that electrons with momenta exceeding the 'free electron' limit of ra e ca§/2 can be 
produced when a laser pulse and a longitudinal electric field interact with an electron via a non- 
wakefield mechanism. The mechanism consists of two stages: the reduction of the electron dephasing 
rate 7— p x /m e c by an accelerating region of electric field and electron acceleration by the laser via the 
Lorentz force. This mechanism can, in principle, produce electrons that have longtudinal momenta 
that is a significant multiple of m e ca^/2. 2D PIC simulations of a relatively simple laser-plasma 
interaction indicate that the generation of super-ponderomotive electrons is strongly affected by this 
'anti-dephasing' mechanism. 



The production of highly energetic electrons in the in- 
teraction of ultra-intense laser pulses with plasmas pQ is 
an essential feature of laser-plasma physics that under- 
pins a wide variety of topics ranging from wakefield accel- 
eration [2] through to Fast Ignition inertial confinement 
fusion [3 . For certain topics the production of electrons 
with the highest energies possible is a matter of specific 
interest. Examples include laser-driven ion acceleration 
schemes based on energetic electrons [Jj (e.g. Target Nor- 
mal Sheath Acceleration and those closely related mech- 
anisms), x-ray generation [5], strong field physics [6], and 
positron production in high Z targets [7]. 

The question as to how to reach high energies, and in 
particular how to exceed what might be termed the 'free 
electron' forward momentum limit of m e ca^/2 (8[ [9] is 
therefore one of general interest. The wakefield scheme 
is one such route to producing high energies, although 
does not directly involve the laser field, and is not very 
well suited to producing high currents of energetic elec- 
trons (c.f. laser interactions with dense plasmas [T0HT4] ) . 
Mechanisms that can breach the 'free electron' limit and 
which directly involve the laser field seem to be more 
subtle, such as the 'Direct Laser Acceleration' scheme 
that takes place in the ion channel produced by trans- 
verse ponderomotive expulsion of electrons by the laser 
pulse. This was first analyzed by Pukhov et al. [15 , and 
was more recently re-analyzed by Arefiev et al. [16]. In 
what follows we define the term 'super-ponderomotive' 
electron to mean an electron with forward momentum 
exceeding m e ca^j2 

In Ref. [16] , a general mechanism of producing super- 
ponderomotive electrons was put forward. The interac- 
tion of an electron with a laser field in vacuum will have 
an integral of motion of the form 7 — p x /m e c = R |8] [9]. 
It is also the case that 7 — p x /m e c is the dephasing rate 
of the electron, and the 'free electron' momentum limit 
arises from R = 1 in the absence of any fields apart 
from the laser field. Arefiev et al. showed that, in the 
specific case of the ion channel, the transverse electro- 



static field in the ion channel can reduce the dephasing 
rate, and thus super-ponderomotive electrons can be pro- 
duced. However there is no a priori reason why only a 
transverse electric field can reduce R below unity. 

In this Letter we show that this very general mech- 
anism extends to the longitudinal electric field as well, 
i.e. collective electric field in the direction of laser prop- 
agation. Longitudinal electric fields are also naturally 
established in laser-plasma interactions by ponderomo- 
tive displacement of electrons, so they are a clear 'next 
candidate' for extra fields that could reduce the dephas- 
ing rate. We show that super-ponderomotive energies 
can result from electron interactions with spikes of rela- 
tively weak longitudinal electrostatic field. In contrast 
with the wakefield acceleration, the axial acceleration 
in this case is insignificant in terms of the direct en- 
ergy gain. Instead, the role of the longitudinal field is 
a reduction of dephasing that subsequently allows for an 
extended interaction with the laser field and leads to a 
significant energy enhancement. The strongest increase 
in energy occurs if the interaction with the longitudi- 
nal field is terminated as the electron passes through a 
zero in the vector potential of the laser field. This in- 
dicates that there are certain parallels with other rela- 
tivists laser-plasma phenomena where vector potential 
zeroes are critical [T7HT9] . This mechanism can work in 
a region of underdense to near-critical plasma formed in 
front of a dense target. Therefore, it may be the case that 
the mechanism is partly responsible for the production of 
highly energetic electrons in current and extant experi- 
ments. To show that this mechanism can naturally occur 
in laser-target interactions, we also present a 2D particle- 
in-cell (PIC) simulation of an extended underdense tar- 
get irradiated by an intense laser beam. We demonstrate 
that a self-consistently generated quasi-static longitudi- 
nal electric field at the edge of the plasma significantly 
reduces the dephasing rate and enables production of 
super-ponderomotive electrons by the laser beam. 

Consider the dynamics of a single electron in an essen- 
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tially ID configuration in which it interacts with a plane 
electromagnetic wave described by the vector potential, 



A = [0, 0, A] = [0, 0, A cos (lj l t) 



(1) 



where r = t — x/c, and ujl is the frequency of the field. 
The electric and magnetic fields are related to the vector 
potential via E = —d t A and B = V x A, so the electric 
field of this wave is polarized in the z-direction. We also 
consider the case where a longitudinal electric field, E x , 
is present. The equations of motion of the electron that 
need to be considered are: 



dp x 
dt 



-eE x + ev z B y , 
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From the definition of r (r 
entiate to obtain, 



t — x/c), one can differ- 
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and this can then be used to write the field components 
as E z = — d T A, By = (l/c)d T A. These can then be 
used to obtain, p z = eA, from Eq. [4j which is one of the 
key integrals of motion. In the absence of E x , another 
integral of motion is obtained from Eq.s[2]and[5j namely 
7 — p x /m e c = 1 (assuming that the electron is initially 
at rest). Using this, one obtains p x = e 2 A 2 /2m e c in the 
E x = case (i.e. the 'free electron' case). If, however, 
E x = — E (where E is a positive constant over some 
region) then we instead have, 
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and from this we can see that 7 — p x /m e c < 1. We can 
now re- write Eq. [2] as, 



where 
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dt 
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(8) 



(9) 



Prom Eqs. [8] and [9] one can see that the effect of the 
accelerating electric field will not only be direct acceler- 
ation of the electron (similar to wakefield acceleration), 



but it will also be a reduction of the dephasing rate R. 
As a result, the electron will gain additional energy from 
the laser field above that obtained in the free electron 
case, i.e. it can produce super-ponderomotive electrons. 
Equation [5] also emphasizes that the 'j x B' force is not 
entirely separated from the longitudinal electric force, as 
the two are linked through the dephasing rate. 

If the electric field has a limited spatial extent, then 
after passing through this spike one will have E x = 0, 
but it will still be the case that 7 — p x /m e c = R, and 
Eq. [8] can then be directly integrated to give, 



Px 



Px 



m e c cr 



,*2 



R 



(10) 



where we have normalized the vector potential via a = 
eA/m e c, and where p* and a* are the longitudinal mo- 
mentum and normalized vector potential immediately af- 
ter the interaction with the spike in the electric field. It 
is clear that the largest effect will be obtained if the de- 
phasing rate R is significantly reduced and the region of 
interaction with E x terminates close to a zero in the vec- 
tor potential (a* < ao). At a* = 0, we have p z = and 
we immediately find that the reduced dephasing rate in 
this case is given by 



R 



7 - pl/m e c = \J (p*/m e c) 2 + 1 - pl/m e c. (11) 



To significantly decrease the dephasing rate (R < 0.5), 
the longitudinal momentum following the interaction has 
to be relativistic. Specifically, assuming that p x /m e c ^> 
1, we find directly from Eqs. 



R: 



and 



PX ~P; 



([lip and ([10]) that 

m e c/ 2p% 

x (l+« 2 ). 



(12) 



(13) 



Therefore, the axial momentum can be enhanced by as 
much as a factor of 2p*/m e c compared to the 'free elec- 
tron' limit. Equation (13) confirms that the longitu- 



dinal momentum can be significantly enhanced even if 
the change in the longitudinal momentum during the in- 
teraction with the static field is relatively small (a§ ^> 
p*/m e c> 1). 

To estimate the amplitude of the static field required to 
reduce the dephasing rate, R = 7 — p x /m e c, well below 
1, we again use the assumption of a region of uniform 
electric field. Making use of Eq.s[6]and[7]we have dr /dt = 
R/7. From this the two key equations to consider are 
Eqj|and, 



dR 

~ctt 



eE 



R 



m e c R 



(14) 



From Eq{l4] we can see that at the zeros in the vec- 
tor potential (p x ~ 0) we can achieve a rapid fall in 
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R, i.e. dR/dt = —eE/m e c. At the peaks of the laser 
field (p x « e 2 A 2 /2m e c) the reduction in dephasing will 
be much slower, i.e. dR/dt « —(eE/m e c)(2R/a 2 ) ). One 
therefore expects, in general, for the largest drops in de- 
phasing to occur around zeroes in the vector potential. 
To significantly decrease the dephasing rate (R < 0.5) in 
a single spike around a zero in the vector potential, one 
would estimate St « 1/(clo(jJl) which implies E « El/2tt, 
where El is the amplitude of the lase field. With a 
more extended electric field, the actual field magnitude 
required will be significantly less. 

The insights gained from this analysis can be verified 
by direct numerical integration of Eq.s [2}|4] along with 
dx/dt = p x /jm e for particular laser fields and choices of 
E x . Here we show the result of a calculation performed 
with a laser field defined (in terms of normalized vector 
potential) by, 



clq cos(co>lt) exp 



-(x — ct — xo) 2 



(15) 



where ao = 10, A =l/im, xo = 6ctL, and tL =40 fs. 
The electron is initially at rest at the origin. A constant 
longitudinal electric field had a magnitude equal to E x = 
—0.1E L over the region 142/im> x > 147/im. The E x 
component is zero at all other points. The change in the 
electrostatic potential across this region is approximately 
16 MeV, which would result in accelerating the electron 
to p x /m e c ^31 without the laser field. 

The results of this calculation shown in Fig. [T] demon- 
strate that the effect of the longitudinal field is to drive 
the electron onto a strongly super-ponderomotive trajec- 
tory. It is evident from the plot of a at the electron loca- 
tion that the interaction with the longitudinal field lasts 
less than a single oscillation of the laser field (the red 
segment of the curve). The immediate effect of the axial 
electric field in this specific case is negative, as the ax- 
ial momentum decreases during the interaction compared 
to the axial momentum calculated with E x = 0. How- 
ever, the interaction also leads to a considerable drop in 
7 —p x /m e c and it terminates close to a zero in the vector 
potential. As a result, a subsequent interaction with the 
laser field leads to a significant longitudinal acceleration, 
with the peak momentum in the excess of 1000 m e c. This 
result is consistent with Eq. ( [To] ) since p*/m e c ~ 27.52, 
R « 0.042, a = 10, and a*/a « 0.116. 

Note that if the longitudinal electric field does not ac- 
celerate the electrons in the direction of the laser pulse, 
then the effect of the longitudinal electric field will be 
instead to reduce the electron momentum. In actual 
laser-plasma interactions, there is a considerable degree 
of complexity, and this means that a full assessment of 
this anti-dephasing mechanism requires a fuller model for 
a more complete assessment. 

Nonetheless the insights gained from studying single 
electron motion are of great benefit when the results of 
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FIG. 1. Results of the case A calculation for single electron 
motion (red line). Region of accelerating E x located over 
142 > x > 147/im. Plot shows (a) longitudinal component 
of electron momentum, (b) vector potential at the electron 
location; (c) the dephasing rate, 7 — p x /m e c, and (d) the 
relativistic 7- factor. Black line in upper panel is the same 
plot for the case without E x . 



2D Particle-In-Cell (PIC) simulations are analyzed in de- 
tail. A 2D PIC simulation was carried out using the 
Plasma Simulation Code (PSC), in which a laser pulse 
was normally incident significantly underdense hydro- 
genic plasma at a density of 8xl0 25 m -3 . The length 
and width of the slab being 200/im and 160/im respec- 
tively, and the simulation domain was 300/im by 200/im. 
The laser pulse had an ao of 10, a FWHM width of 8/im, 
and a pulse duration of 500 fs. Denoting the two axes 
of the simulation domain as x (laser direction) and y, 
the electric field of the laser pulse was polarized in the z 
direction. 

This pulse is significantly longer than the characteris- 
tic time of electron response in the plasma, so that the 
laser eventually creates a quasistatic channel shown in 
Fig. [2j The channel has a coaxial structure with a posi- 
tively charged center and negatively charged walls. Dur- 
ing the formation of the channel, the electrons from in- 
side the channel are expelled radially by ponderomotive 
pressure, which causes charge separation. The resulting 
radial electric field counteracts the expelling force, allow- 
ing the electrons to remain in an equilibrium bunched 
on the perephery of the laser beam. Such channels and 
the corresponding transverse electric field are routinely 
observed in simulations of laser interactions with under- 
dense plasmas [20-23 . However, the fact that such a 
coaxial structure would also produce a quasistatic axial 
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electric field at the channel entrance (see Fig. [2]) has been 
underappreciated. 

In order to examine the role played by the axial field, 
we have performed a search on the electron particle data 
for super-ponderomotive electrons for which the dephas- 
ing rate satisfied 7 —p x /m e c <0.05 and x <50/im at least 
at some point during the electron trajectory. In Fig. [2] 
we show a trajectory, axial momentum, and dephasing 
for just one such electron. Near the front of the slab, 
there is a region with a strong quasi-static negative E x . 
Clearly, the electron interaction with this field (shown on 
all plots with a red segment) launches the electron onto a 
super-ponderomotive trajectory (the subsequent acceler- 
ation is shown with a light-blue curve). There is virtually 
no self- focusing of the laser pulse in this region {clq ~10), 
so the free-electron limit for 7 is 50. The electron how- 
ever achieves a peak 7 of nearly 150, which exceeds the 
free-electron limit by a factor of three. The accelera- 
tion is preceded by a massive drop in the dephasing rate 
that occurs during the interaction with E x when the elec- 
tron moving against the laser beam is turned around and 
pushed forward. The time evolution of the axial momen- 
tum further emphasizes that this is a two-stage process, 
since no significant axial acceleration occurs directly dur- 
ing the interaction with the axial field. 

It must be pointed out that the electric field of the 
laser is polarized out of the plane of the simulation do- 
main, which eliminates forced driving of betatron oscil- 
lations [15] as a possible explanation for the observed 
significant energy gain. Moreover, we observe no ampli- 
fication of the transverse oscillations across the channel, 
which indicates that the observed effect is not related to 
the parametric amplification of betatron oscillations [16 . 
Figure [2] clearly shows that this event is quite prompt, so 
the underlying mechanism must be able to produce the 
observed behavior without any gradual build-up. The re- 
duction of dephasing by acceleration in the longitudinal 
electric field clearly satisfies this key criterion. 

Later on it can be seen that there is a decline in the 
electron momentum. This illustrates the concern, stated 
earlier, that the collective fields can also act to increase 
the dephasing rate. Note that in this simulation there 
is no dense foil which will interrupt the acceleration pro- 
cess. We have therefore made no attempt to determine 
the optimal conditions for exploiting this anti-dephasing 
mechanism, and this is matter for future work as these 
conditions will be highly dependent on the specifics of 
the laser and target parameters. 

In this Letter we have shown that, according to the 
single electron equations of motion, electrons can be ac- 
celerated to momenta in excess of m e ca^j2 by an 'anti- 
dephasing' mechanism in which a brief acceleration by a 
longitudinal electric field that is simulataneously present 
with the laser pulse reduces 7 — p x /m e c significantly 
below unity. This was verified by direct numerical in- 
tegration of the equations of motion and then the ef- 
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FIG. 2. Snapshots of time- averaged axial electric field, elec- 
tron density, and time evolution of the electron axial momen- 
tum and the dephasing from a 2D PIC simulation. 



feet was confirmed in a fully self-consistent simulation 
of laser-plasma interaction (laser pulse propagating in 
under-dense plasma). This therefore shows that the anti- 
dephasing mechanism both exists and can play a very im- 
portant role in the production of super-ponderomotive 
electrons in laser-plasma interactions. This mechanism 
is complimentary to the mechanism of the parametric 
amplification of betatron oscillations, so that the combi- 
nation of the two can produce super-ponderomotive elec- 
trons with energies exceeding what is predicted here and 
in Ref. [16]. The work presented here also shows that 
one cannot simply split electron motion into an indepen- 
dent 'wakefield-like' component and a 'free electron like' 
component, when in fact the anti-dephasing mechanism 
is due to a subtle interaction between the two. 
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